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Abstract. We consider the linear dissipative Boltzmann equation describing inelastic 
interactions of particles with a fixed background. For the simplified model of Maxwell 
molecules first, we give a complete spectral analysis, and deduce from it the optimal rate 
of exponential convergence to equilibrium. Moreover we show the convergence to the heat 
equation in the diffusive limit and compute explicitely the diffusivity. Then for the physical 
model of hard spheres we use a suitable entropy functional for which we prove explicit 
inequality between the relative entropy and the production of entropy to get exponential 
convergence to equilibrium with explicit rate. The proof is based on inequalities between 
the entropy production functional for hard spheres and Maxwell molecules. Mathematical 
proof of the convergence to some heat equation in the diffusive limit is also given. From 
the last two points we deduce the first explicit estimates on the diffusive coefficient in the 
Fick's law for (inelastic hard-spheres) dissipative gases. 

1. Introduction. The linear Bohzmann equation for granular particles models the dy- 
namics of dilute particles (test particles with negligible mutual interactions) immersed in a 
fluid at thermal equilibrium that undergo inelastic collisions characterized by the fact that 
the total kinetic energy of the system is dissipated during collision. Such an equation in- 
troduced in 1 16, 21 , 15 1 provides an efficient description of the dynamics of a mixture of 
impurities in a gas ||T3llT0l . Assuming the fluid at thermal equilibrium and neglecting the 
mutual interactions of the particles, the evolution of the distribution of the particles phase 
is modelled by the linear Boltzmann equation which reads 

dtf + vVJ = Qif), (1.1) 

with suitable initial condition f(x,v,Q) = fo{x,v), (x,v) e X R^. Here above, the 
collision operator Q{-) is a linear scattering operator given by 

Q(/) = S(/,Mi), 
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where !B{-, ■) denotes the usual quadratic Boltzmann collision operator for granular gases 
(cf. 191 for instance) and Aii stands for the distribution function of the host fluid which is 
assumed to be a given Maxwellian with bulk velocity Mi and temperature ©i (see Section 2 
for details). Notice that we shall deal in this paper with the collision operator Q correspond- 
ing to hard-spheres interactions as well as with the one associated to Maxwell molecules 
Interactions . The inelasticity is modeled by a constant normal restitution coefficient. The 
main goals and results of this paper are the following: 

(1) First, we explicit the exponential rate of convergence towards equilibrium for the so- 
lution to the space-homogeneous version of (II. lb (for both Maxwellian molecules and 
hard-spheres interactions) through a quantitative estimate of the spectral gap of Q. 
It is computed in an explicit way for Maxwell molecules, and estimated in an explicit 
way for hard-spheres, see Theorem 13 . 8 1 (to gether with its Corollarv l3.9| for its conse- 
quence on the asymptotic behavior of space-homogeneous solutions), where /.(max is 
defined in Theorem 13. 21 and the constant C* is detailed in Remark [33] 

(2) Second, we investigate the problem of the diffusion approximation for (II. lb . Pre- 
cisely, we show that the macroscopic limit g of Eq. (II. lb in the diffusive scaling 
is the solution to some (parabolic) heat equation (see Proposition 14.31 together with 
Theorems 14.41 and 14. 9b . When dealing with Maxwell molecules interactions, the 
diffusivity of this heat equation can then be computed explicitly (see Theorem 14.41 
together with the computation of Remark |43] for the diffusivity). This is no more 
the case for the equation corresponding to hard-spheres interactions but we provide 
some new quantitative estimates on it (see Theorem |4.9| together with the estimate of 
Proposition l4.10b . 

Concerning point (1), it is known from ilT6l 1211 [TSl that the linear collision operator 
admits a unique steady state given by a (normalized) Maxwellian distribution function M 
with bulk velocity U\ and temperature 0* < ©i. Moreover, thanks to the spectral analysis 
of Q performed in 1 1 1 (in the hard-spheres case), the solution to the space-homogeneous 
version of ( II. lb is known to converge exponentially (in some pertinent norm) towards 
this equilibrium M as times goes to infinity. This exponential convergence result is based 
upon the existence of a positive spectral gap for the collision operator Q and relies on 
compactness arguments, via Weyl's Theorem. Because of this non constructive approach, 
at least for hard-spheres interactions, no explicit estimate on the relaxation rate were 
available by now. It is one of the objectives of this paper to fill this blank. It is well-known 
that the kinetic description of gases through the Boltzmann equation is relevant only on 
some suitable time scale [9", 1 1 1. Providing explicit estimates of the relaxation rate is the 
only way to make sure that the time scale for the equilibration process is smaller than the 
one on which the kinetic modeling is relevant. Another motivation to look for an explicit 
relaxation rate relies more on methodological aspects. Compactness methods do not rely 
on any physical argument and it seems to us more natural to look for a method which relies 
as much as possible on physical mechanisms, e.g. dissipation of entropy. 

Precisely, the strategy we adopt to treat the above point (1) is based upon an explicit 
estimate of the spectral gap of Q. For Maxwell molecules interactions, we use the Fourier- 
based approach introduced by Bobylev ||6l for the study of the linearized (elastic) Boltz- 
mann equation and then used for the study of the spectrum of the linearized inelastic col- 
lision operator in [7 1, and we provide an explicit description of the whole spectrum of this 
linear scattering collision operator . Then, to treat the case of hard-spheres interactions, 
our method is based upon the entropy-entropy-production method. Precisely, we show that 
the entropy production functional (naturally associated to the L^(A1~^) norm) correspond- 
ing to the hard-spheres model can be bounded from below (up to some explicit constant) 
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by the one associated to the Maxwell molecules model (Proposition [33]l. Note that such 
a comparison between entropy production rates for hard-spheres and Maxwell molecules 
interactions is inspired by the approach of Q which deals with the linearized (elastic) 
Boltzmann equation. In the present case, the method of proof is different and simpler, be- 
ing based upon the careful study of a convolution integral. Such an entropy production 
estimate allows us to prove, via a suitable coercivity estimate of Q (Theorem l3.8l l, that any 
space-homogenous solution to ( II. Il l converges exponentially towards equilibrium with an 
explicit rate that depends on the model under investigation. 

Concerning now point (2), various attempts to derive hydrodynamic equations from the 
dissipative nonlinear Boltzmann equation exist in the literature, mostly based upon suitable 
moment closure methods ||23] H] |5l or on the study of the linearized version of the Boltz- 
mann equation around self-similar solutions (homogeneous cooling state) fT, '8^ in some 
weak inelastic regime. Dealing with the linear Boltzmann equation (II. Il l, hydrodynamic 
models describing the evolution of the momentum and the temperature of the gas have been 
obtained in 1 10] as a closed set of dissipative Euler equations for some pseudo-Maxwellian 
approximation of Q. Similar results have been obtained in 1 13 1 where numerical methods 
are proposed for the resolution of both the kinetic and hydrodynamic models. The work Q 
proposes two closure methods, based upon a maximum entropy principle, of the moment 
equations for the density, macroscopic velocity and temperature. These closure methods 
lead to a single diffusion equation for the hydrodynamical variable. In the present paper, we 
shall discuss the diffusion approximation of the linear Boltzmann equation dl.lb with the 
main objective of providing a rigorous derivation of the Pick's law for dissipative gases 
and an estimate on the diffusive coefficient. Recall that the diffusion approximation for 
the linear Boltzmann equation consists in looking for the limit, as the small parameter e 
goes to 0, of the solution to the following re-scaled kinetic equation: 

edtMt,x,v) + V ■ VJS,x,v) = ^Qif,){t,x,v), (1.2) 

with suitable initial condition. We consider indeed here the Navier-Stokes scaling, namely, 
we assume the mean free path to be a small parameter A = e ^ 1 and, at the same time, 
we rescale time as t —> t/e in order to see emerging the diffusive hydrodynamical regime 
(and not the Euler hydrodynamical description, which would be a trivial transport equation 
in our case). Performing a formal Hilbert asymptotic expansion of the solution allows us 
to expect the solution to converges towards a limit / with Q{f) - 0. Therefore, the 
expected limit of is of the form f{t, x, v) = g{t, x)Ai{v), and the diffusion approximation 
problem consists in expressing Q{t, x) as the solution to some suitable diffusion equation. 
Actually, standard approach consists in using the continuity relation 

dtQ{t,x) + div.t/(t, x) = 0, 

between the density g and the current vector i{t, x) together with a suitable Pick's law that 
links the current / to the gradient of g: 

i{t,x) = -DW-,g{t,x) 

for some suitable diffusion coefficient (diffusivity) D > which depends on the kind of 
interactions we are dealing with. For Maxwell molecules interactions, the expression of 
the diffusivity can be made explicit while this is no more the case when dealing with hard- 
spheres model. The method we adopt for the proof of the diffusive limit follows very 
closely the work of P. Degond, T. Goudon and F. Poupaud lfT2l . Though more general than 
ours since it deals with models without detailed balance law, the study of [121 is restricted 
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to the case of a collision operator for which the collision frequency is controlled from above 
in a way that excludes the case of physical hard-spheres interactions (recall that, for hard- 
spheres, the collision frequency behaves asymptotically like (1 + \v\) 11]). Actually, the 
analysis of fT2l can be make valid under the only hypothesis that the coercivity estimate 
obtained in Theorem 13.81 (and strenghten in Theorem 13. 10b holds true. Precisely, such a 
coercivity estimate of Q allows to obtain satisfactory a priori bounds for the solution to the 
re-scaled equation ( 11.21 ). We are then able to prove the weak convergence of the density 
and current of the solution towards suitable limit density g and limit current j. It is 
also possible, via compensated-compactness arguments from lfT4l [iTl l22l to prove strong 
convergence result in norm. 

The organization of the paper is as follows. In Section|2]we present the models we shall 
deal with as well as some related known results we shall use in the sequel. In Section[3]we 
perform the computations of the spectrum in the Maxwell molecules case and then prove 
the crucial entropy production estimates in the hard-spheres case (from which we deduces 
the explicit convergence rate to equilibrium for the space-homogenous version of ( 11.11 )). 
Section|4]is dealing with the above point (2). We first prove a priori estimates valid for both 
models of hard-spheres and Maxwell molecules and based upon the coercivity estimates 
obtained in Section [3] Then, we deal separately with the cases of Maxwell molecules 
and hard-spheres proving for both models the convergence towards suitable macroscopic 
equations, providing for Maxwell molecules an explicit expression of the diffusivity, and 
for hard-spheres explicit estimates on it. 

2. Preliminaries. 

2.1. The model. As explained in the introduction, we shall deal with a linear scattering 
operator Q given by 

Q{f) = ^ f B{q,n)[}Bf{x,v^,t)Mi{zv^)-f{x,v,t)Mi{w)]dwdn (2.1) 

where A is the mean free path, q - v - w is the relative velocity, and if* are the pre- 
coUisional velocities which result, respectively, in v and w after collision. The main feature 
of the binary dissipative collisions is that part of the normal relative velocity is lost in the 
interaction, so that 

(c* -IV*) -n = -e{v - w) ■ n, (2.2) 

where n e is the unit vector in the direction of impact and < e < 1 is the so-called 
normal restitution coefficient. Generally, such a coefficient should depend on (v, w) but, 
for simplicity, we shall only deal with a constant normal restitution coefficient e. The 
collision kernel B{q, n) depends on the microscopic interaction (see below) while the term 
/b corresponds to the product of the Jacobian of the transformation {Vi,, w*) {v, w) with 
the ratio of the lengths of the collision cylinders |9|. Note that in such a scattering model, 
the microscopic masses of the dilute particles m and that of the host particles nii can be 
different. We will assume throughout this paper that the distribution function M\ of the 
host fluid is a given normalized Maxwellian function: 

where U\ e is the given bulk velocity and ©i > is the given effective temperature 
of the host fluid. For particles of masses m colliding inelastically with particles of mass 
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mi, the restitution coefficient being constant, the expression of the pre-colHsional velocities 
{Vi,, w*) are given by 19 , 211 

1-B 1-6 
v^=v-la ^_^^ {q-n)n, = w + 2{1 - a) ^ _^^ {q ■ n) n; (2.4) 

where q = v - w, ais, the mass ratio and jS denotes the inelasticity parameter 

e(0,l), iS = i-^e [0,1/2). 



ni + nii 



We shall investigate in this paper several collision operators corresponding to various inter- 
actions collision kernels. Namely, we will deal with 

• the hnear Boltzmann operator for hard-sphere interactions Q - Qhs for which 

1 

B{q, n) = Bhs((?, n)^\q- n\, and /g,^ /hs = ^; 

• the scattering operator Q - Qmax corresponding to the Maxwell molecules approxi- 
mation for which 

^ 1 |l7 — w\ 

B(q,n) = Bmax{q,n) = \q-n\, q = q/\q\, and Jb^^ = /ma 



\v* - w* 



It will be sometimes convenient to express the collision operator Q in the following weak 
form: 

<ip,Q{f)>=-^ f B{q,n)f{v)Mi{w)[ip{v*)-4>iv)]dvdivdn (2.5) 

ZnA Jr3x]r3xS2 

for any regular ip, where (v*, zv*) denote the post-coUisional velocities given by 

V* ^v-2a{l-^){q-n)n, w* ^ iv + 2{1 - a){l - ^){q ■ n)n. (2.6) 

In particular, one sees that the dissipative feature of microscopic collision is measured, at 
the macroscopic level, only through the parameter: 

K = a(l-|3) = -(l + e)e(0,l) 

appearing in the expression of v*. Accordingly, the macroscopic properties of Q are those 
of the classical linear Boltzmann gas whenever k = 1/2 (which is equivalent to v* = v). 
It can be shown for both cases that the number density of the dilute gas is the unique 
conserved macroscopic quantity (as in the elastic case). In contrast with the nonlinear 
Boltzmann equation for granular gases, the temperature, though not conserved, remains 
bounded away from zero, which prevents the solution to the linear Boltzmann equation to 
converge towards a Dirac mass. 

Moreover let us remark that from the dual form we see that the collision operator in fact 
depends only on two real parameters m\ and k (plus Ui of course) and not Ui plus three 
parameters a, /3, ni\ as a first guess would suggest. 

2.2. Universal equilibrium and H-Theorem. A very important feature of these inelastic 
scattering models is the existence (and uniqueness) of a universal equilibrium, that is in- 
dependent of range of the microscopic-interactions (that is of the collision kernel B) and 
depending only on the parameters m\, k and U\. Precisely, the background forces the sys- 
tem to adopt a Maxellian steady state (with density equal to 1): 
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Theorem 2.1. The Maxwellian velocity distribution: 

.... I m ( m(v-uif-\ 



with 

. (l-a)(l-iS) 
® = 1 TT-^®^ (2.8) 

;i the unique equilibrium state ofQ with unit mass. 

Note that this universal equilibrium is coherent with the remark that the collision only 
depends on m\ and k (and U\) from the dual form, since 

m _ mi 1 - K 

e* ~ ©l K ■ 

This explicit Maxwellian equilibrium state allows to develop entropy, spectral and hy- 
drodynamical analysis on both the models in the same way. First, from 1 18, 19| and ifTSl . 
the existence and uniqueness of such an equilibrium state allows to establish a linear ver- 
sion of the famous H-Theorem. Precisely, for any convex C^-function O : IR""" M, the 
associated so-called H-functional (relatively to the equilibrium Ai) 



H«(/IM) 



is decreasing along the flow of the equation ( 11.11 ) (this is the opposite of a physical en- 
tropy), with its associated dissipation functional vanishing only when / is co-linear to the 
equilibrium Ai: 

Theorem 2.2 (Formal H-Theorem). Let f{t, v) ^ be a space homogeneous solution 
then we have formally 

^H^if{t)\M) = ^^^f^^^'^'^'^'l^) ^ {t> 0). (2.10) 

The application of the H-Theorem with <I'(x) = {x - 1)^ suggests the following Hilbert 
space setting: the unknown distribution / has to belong to the weighted Hilbert space 
L2(A1-1) = L'^(Rl; M~^{v)dv). Consequently, one defines the Maxwell molecules and 
hard spheres collision operators, associated to the mean-free path A - 1, with their suitable 
domains, as follows: 

l^(Xhs) c L2(M-1), Range(Xhs) c L^M-'), 
jXhs/ = Qhsf for any / e ^(Xhs), 



and 



Precisely, 



( ^(£max) = LHM-^), Range(Xmax) c L2(M-1), 
I -Cmax / = Qmax / for any / € ^(Xmax)- 



where fftis is the collision frequency associated to the hard-spheres collision kernel: 

Ohs{^) = J- f 1^ ■ n\Mi{zv) dzvdn, v e R^. 
Note that crhs is unbounded [|1||: there exist positive constants Vq, Vi such that 

Vo(l + \V- Mil) < 0'n%{v) < Vi(l + \V- Mil), G 



INELASTIC SCATTERING BOLTZMANN MODELS 



7 



For this reason, ^(Xhs) L?-{Ai ^). On the contrary, the collision frequency cJmax associ- 
ated to the Maxwell molecules collision kernel. 



Omaxi'v) = I?' n\Mi{zv)dwdn 



is independent of the velocity v and Xmax is a bounded operator in L^{Al~^). We recall (see 
ifTl ) that Xhs is a negative self-adjoint operator of L^{M~^). Moreover, let us introduce the 
dissipation entropy functionals associated to Xmax and Xhs- 

l)max(/)-- r £ma.{f){v)f{v)M-\v)dv, f e LHM-') 

and 

£»hs(/):=- r Xhs(/)(t')/(i')M-Hi^)dz;, /ei^(Xhs). 

Note that, by virtue of ( 12.101 ). if /(f) denotes the (unique) solution to ( II. lb in L^(yV(~^) for 
hard-spheres interactions, then, with the choice 0(x) - {x - 1)^, 

^11/(0 - = ^H<p(/(OIM) = -2 2)hs(/(0)- (2.11) 

The same occurs for £)max(/(f))- This is the reason why we are looking for a control 
estimate for both D^sif) and Dmaxif) with respect to the L?-{Ai~^) norm of /. It will be 
useful to derive an alternative expression for both ©max and Dhs- 



Proposition 2.3. For any f e ^(Xhs), 

f(v*) fiv) f 

■'^ ' ' Mi{w)M(v) dw dv dn > 0. 



yv((z;*) yv((i;) 



In the same way, 

f{v*) f{v) 



forany f eh^iMr^). 



M{v*) M{v) 



M\{w)M{v) dw dv dn ^ 0. 



Proof. The proof is a straightforward particular case of the above H-Theorem. Precisely, 
let us fix / e ^(Xhs). and set / = gM then one has 

r £Uf)fM-^ dv= { ^^M{v)Mi{w)g{v) {g{v*) - g{v)] dw dv dn. 

Jr3 Jr3xR3xS2 ^™ 

Moreover, for any (p e f^(Xhs)' since M is an equilibrium state of (Shs^ 

r \q-n\ , ^ 

- — M{v)Mi{w) [(p{v*) - (p{v)\ dw dv dn = {£^,s{M), (ph^nAv) = 0- 

Jr3xR3xS2 '^^'^ 

It is easy to deduce then that 
r Xhs(/)/M-'dz;= r ^■^M{v)Mi{w)g{v*)[g{v)-g{v*)] dwdvdn. 

Jr3 Jr3xR3xS2 ^™ 

Finally taking the mean of these two quantities leads to the desired result. The same rea- 
soning also holds for ©max- n 
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3. Quantitative estimates of the spectral gap. In this section, we strengthen the above 
resuh in providing a quantitative lower bound for both £)max(/) and Dhs(/)- The estimate 
for Dmaxif) is related to the spectral properties of the collision operator Xmax while that 
of £)hs(/) relies on a suitable comparison with Dmaxif)- From now on, (•, •) denotes the 
scalar product of L^{M~^) and S(X), shall denote respectively the spectrum and 

the point spectrum of a given (non necessarily bounded) operator in L^(A1~^). 

3.1. Spectral study for Maxwell molecules. We already saw that the operator Xmax : 
L^{M~^) L^{M~^) is bounded, and it is easily seen that Xmax splits as 

-Cmax/ = -Cmax/ - /(^) 
where Xmax is compact and self-adjoint (the proof can be done similarly as in [1 J) and we 
used that the collision frequency ffmax associated to the Maxwell molecules is constant. 
Moreover, since 

<Xmax/,/)<0, V/eL2(M-i), 
the operator is easily seen to generate a C°-semigroup of contractions, and it is known that 
S(Xmax) C (-00,0]. Finally, since Id + Xmax ~ -£max is ^ positive, self-adjoint compact 
operator, one sees that the spectrum of Xmax is made of a discrete set of eigenvalues with 
finite algebraic multiplicities plus possibly {-1} in the essential spectrum, with 

Sp(Xmax)c(-l,0] (3.1) 

and where the only possible accumulation point is {-1}. Clearly, since Xmax(AI) - 0, 
Ao,o := is an eigenvalue of Xmax with eigenspace given by Span(AI). There are other 
eigenvalues of Xmax of peculiar interest. Namely, for any / e L^{Ai~^), the weak formula- 
tion ( 12.5b yields 

I {v- ui)£max{f) dv = ^ ( Iq- n\f{v)Mi{w) [v* - v] dvdwdn 

I f{v)dv I Miiw)div I \q- n\{q ■ n)ndn. 

Jr3 Jr3 Js2 

Using the fact that I \q ■ n\{q ■ n)ndn = nq, one has 

I {v - ui)£rnax{f) dv = -a{l - I f{v)dv I {v - w)Mi{w) dw 

= -a(l-jS) ( {v-ui)f{v)dv. 
Jr? 

The operator Xmax being self-adjoint in L^{Ai~^), one obtains the identity 

</,i;max((i^, - uu)M)) = -ail-p){f,{v, - uu)M), i = 1,2,3, V/ e L^M.-'). 

If we denote 

Ao,i = a(l-iS) = Ke(0,l), 
this means that -Ao,i is an eigenvalue of Xmax associated to the momentum eigenvectors 
{Vi — U\j)Ai{v) for any i = 1,2,3. In the same way, technical calculations show that, for 
any/eL2(M-i), 



/R^xR^xS^ 

a(l-jS) 



<Xmax(/),|l'-Mll'M)= f Xmax(/)l^'-Mll'dz; 

Jr3 



= -2a(l - /3)(1 - a(l - )?))(/, \v - mfM) + - ^)'< f, M), 
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and, since Xmax is self-adjoint, setting 

Ai,o = 2a(l - - a{l - |S)) = 2k(1 - k), 
one has {£max{\'^-Ui\^M), f) = -^-ifiilv-Uif-M, f), for any /± span(Al). Equivalently, 

Xmax (£) ~ -Ai,o£, 

where S{v) = (\v - Mip - ^) A1, is the energy eigenfunction, associated to -Aj^o- To 
summarize, we obtained three particular eigenvalues Ao,0/ Ao,i and Ai,o of Xmax associated 
respectively to the equilibrium, momentum and energy eigenfunctions. 

To provide a full picture of the spectrum of Xmax, we adopt the strategy of Bobylev 
||6l based on the application of the Fourier transform to the elastic Boltzmann equation. 
Namely, we are looking for A > such that the equation 

Xmax/ = -A/, / e l2(M-1), / ^ 0, (3.3) 

admits a solution. Applying the Fourier transform J? to the both sides of the above equa- 
tion, we are lead to: 

where / = ^{f). One deduces immediately from the calculations performed in 11271 . that 

^ [Xmax/] (?) = ^ £^ I?- "I [TiC)Mi{n - MiMm] dn 

with ^ = 4/|<£| and 

C = 2a(l - ^){E, ■ n)n, T = 2a(l - /S)(5 • n)n 

while M\ is the Fourier transform of the background Maxwellian distribution, given by 



M(5) = expl-iui-E.- 



2mi 



The fundamental property is that even if the equilibrium distribution Ai does not make the 
integrand of the collision operator vanish pointwise (as it is the case in the elastic case), 
surprisingly it still satisfies a pointwise relation in Fourier variables as was noticed in II2TI . 
A simple computation yields 

M(5) = exp|-ZMi-5-^^ 

and thus one checks easily that, M{£,^) Mi{£,~) = M{£,), for any n e S^, and any £, e R^. 
Thus, following the method of Bobylev [6J, we rescale - (p(5) and define the 

corresponding re-scaled operator ^ 



271 Js2 



Then, Eq. ( I3.3l l amounts to find A > such that the equation ^(p{£,) = -A(p{£,) admits a 
non zero solution q) with 

/ = {<pM-^) e (M"^ ) (3 .4) 

where stands for the inverse Fourier transform. Using, as in the elastic case f6^, p. 
136], the symmetry properties of the operator ^ together with condition (13.4b . one obtains 
that functions of the form 
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Yf „, being a spherical harmonic e N, m = ...,£), are eigenvectors of associated 
to the eigenvalues -A„/ - -A„/{k) where, according to |7|, 

A„ , = 1 - — ^ C s^"^'^' I ^.~.^Vf ] ds (3.5) 
2k{1 - k) Ji_2, ' \ (2 - 2k)s j 

where P/" is the ^-th Legendre polynomial, n e N, ^ e N. Note that, according to ( 13. It . 
A„/ e (0, 1]. Technical calculations prove that the eigenvalues we already found out, 
namely, 

Ao,o = 0, Ao,i = K, and Ai,o = 2k(1 - k), 
do actually correspond to the couples {n, £) = (0, 0) ; (0, 1) and (1, 0) respectively. More- 
over, from the well-known Legendre polynomials property: 

{( + l)Pf+i(x) = {li + l)xPf(x) - i Pf-i(x), X e ]R, ^ > 1, 

one obtains the recurrence formula 

V ^ 2^+1 , e ^ o 

A„,m = — — A„,, + ^^^^^^A„.u-— A„.i,,_i, n,t>Q (3.6) 
where v = 1 - 2k e (-1, 1). Such a recurrence formula together with the relation 

1 1 - v2«+2 



A„,o = 1 - 



n + 1 1 - v2 
allow to prove by induction over ^ e N that 

An+i/ ^ A„/ and A„/+i ^ A„/ for any n e N. 

Consequently, one sees that min{A„/ ; «, /' ^ 0} \ {0} = min{Ai o; Ao,i} which means that 
spectral gap of Xmax is given by 

fmax = min{Ai,o;Ao,i} = min{K;2K(l - k)}. 

Remark 3.1. f^ote that, 

Ao,i < Ai,o <=^ K < 1/2 <;=^ a(l + e)< 1. 

In particular, if mi ^ m then Ao,i ^ Ai^o- Assuming for a while that we are dealing with 
species of gases with same masses ni = m\, then in the true inelastic case (i.e., e < 1) one 
also has Aj^o > A04. This situation is very particular to inelastic scattering and means 
that the cooling process of the temperature happens more rapidly than the forcing of the 
momentum by the background, whereas when K = 1/2 these two processes happen at 
exactly the same speed (Ao,i = Ai,o). Note also that, whenever Ao,i > Ai,o (due to the ratio 
of mass different from 1 ), the smallest eigenvalue corresponds to the momentum relaxation 
and not the energy relaxation anymore. This contrasts very much with the linearized case 
(see Qj. Note also that the first eigenvalues are ordered as illustrated in Fig. 1. 

The above result can be summarized in the following where the last statement follows 
from the fact that Xmax + Id is a self-adjoint compact operator of L^(A1~^). 

Theorem 3.2. The operator — Xmax is a bounded self-adjoint positive operator of L^{M~^) 
whose spectrum is composed of an essential part {+1} plus the following discrete part: 

Sp(-Xmax) = [An/; M, ^ e N) c [0,1) 
where A,,/ is given by ( 13. 51 ). Moreover, Aq^o = is a simple eigenvalue of Ji^ax associated 
to the eigenvector M and — Xmax admits a positive spectral gap 

/"max = inin(Ai,o; Ao,i) = min(K;2K(l - k)}. 
Finally, there exists a Hilbert basis ofL'^(M~^) made of eigenvectors q/— Xmax- 
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Figure 1 . Evolution of Ao,i (black solid line), Ai,o (black dotted line), 
Ai4 (grey dotted line) and Ao,2 (grey solid line) as functions of c e [0, 1] 



3.2. Entropy estimate for Maxwell molecules. The result of the above section allows to 
provide a quantitative version of the H- Theorem. Precisely, for any / e L^(A1~^) orthog- 
onal to At, using the decomposition of both Xmax/ and / on the Hilbert basis of L^(yV(~^) 
made of eigenvectors of -Xmax (see Theorem l3.2l ). it is easily proved that: 

l^maxa) - - J^^ Xmax(/) /M"' Av > fi^ax \\ffn(M-^y V/± Span(M). (3.7) 

It is well-known that such a coercivity estimate allows to obtain an exponential relax- 
ation rate to equiUbrium for the solution to the space homogeneous Boltzmann equation. 
Namely, given /o(z?) e L^{Bp,M~^{v) dv) with unit mass 

r /o(i^)dz; = l, 

let ft be the unique solution of ( II. Il l with initial condition ft=o = /q. According to the 
conservation of mass, it is clear that {ft- Ai) is orthogonal to M (for the L^(]R^, Ai~^ (v) dv) 
scalar product) and, within the entropy language: 

■^H^iftlM) = -2£)max(/0 < -2^maxH<p(/,|M) 
for - (x - 1)^ or equivalently, 

i^j^ (ft-MfM-' dv'j <|j^ (/o-M)2M-i dz;J exp (- ^i^^^t) , Vf ^ 0. 

We obtain in this way an explicit exponential relaxation rate towards equilibrium for the 
solution to the space homogeneous linear Boltzmann equation which is valid for granular 
gases of Maxwell molecules and generalizes a well-known result for classical gases 1 11 1. 
More interesting is the fact that the knowledge of the spectral gap of Xmax allows to recover 
an explicit estimate of the spectral gap of the linear Boltzmann operator for hard-spheres 
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Xhs through a suitable comparison of the entropy production functionals Dhs and Dmax- 
This is the subject of the following section. 



3.3. Entropy estimate for hard-spheres. The goal of this subsection is to show that 
the entropy production functional Dhs for hard-spheres relates to the one for Maxwell 
molecules Dmax- More precisely we shall show that 

Proposition 3.3. The entropy production functionals £)hs <^nd Dmax ore related by: 

DUf) > C* £)max(/), V/ e ^(Xhs) 
for some explicit constant C* depending only on a and jS. 

Remark 3.4. The idea of searching for such an inequality was already present in f2\, 
but here the method of proof is different and simpler: one does not need any triangular 
inequality between collisions, and the proof reduces to a careful study of a convolution 
integral. 

Remark 3.5. Note that in the hard-spheres case, the operator Xhs unbounded. For a 
careful study of its properties (compactness of the non-local part, definition of the associ- 
ated -semigroup of contraction in the Hilbert space L^{M~^)) we refer to \ \\. 

Proof. Let / e i^(Xhs)- We set Ui = in this proof without restriction since this only 
amounts to a space translation. 

We introduce the following parametrization, for fixed n ^ , v = rn + fj, v* = r*n + v, 
w = rzi,n + w, w* - r-^i,*n + w, where r,r*,r-^i, and Vze* are real numbers and v,w are 
orthogonal to n. Simple computations show that 

r„, = — — + 11- ^ 



2a(l - jS) \ 2a(l - jS) 



while 



= I — — - - ^ I r* + r 



^2a(l-/3) a j \a 2a{l-p)) 

Therefore, ri„ and /"a,* only depend on r and r* . Then if we denote 6 the angle between q 
and n, we get from Prop. 12.31 where we set g = 

£'hs(/)=:7- if I \q\ cos e\g(r*n + v) - g{rn-\-v)f 

Mii^wn + w)A\{rn + v) dV dw dr dr* dn 

with 

1^1 = {\v - + (2)c)-2|r* - r\^f 

and 

(27c)-V* - r\ 



vl/2 



COS a = 



where we recall that k = a(l - jS). We split the integral into two parts according to |r* - r| ^ 
fiQ > or \r* - r| < po where Qq is a positive parameter to be determine latter. Using the 
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fact that \q\ > \r* - r\/2K, one has the following estimate for the first part of the integral 
^ [ dn [ dr dr* f \q\ cos 6 M{rn + v) Ml(r^^,n + w) 

^g{r*n + v)- g{rn + u)j dv dw 

f COS yV((rM + v) M\{riun + ■w)^g{r*n + v)- g{rn + u)j dv dw 

which corresponds (up to the multiplicative factor (2k)~^ Qq) to the integral for |r- r*| > j 
corresponding to Maxwell molecules, i.e., 

- I X\\r-r*\>e«]-CUf)fM-^dv>-^ I A-|k-ri>eo)-£max(/)/W("Mi'. (3. 

Concerning now the second part of the integral (corresponding to \r* - r| < Qq), we use 
that \q\ ^ \v -zD\ and we isolate the integration over w: 



8) 



-f A:i|r-r*|<eo)-Chs(/)/A^"'di;^ f ( drdr*\r-r\ 

M{rn + v) Mi{r-,un) ^g{r*n + v)- g{rn + dv 
where we used the fact that, since w is orthogonal to n, 

I nil 

Mi{r-a,n + w) = [-^^^1 Mii'w)Mi{r^n). 

Setting I; - \r* - r\/2K, if one were able to prove that there is a constant C such that 
r \v-zv\Mi{zo) r Mi{w) 

Jr2 {\v - iDp + Qv - zvf + 

uniformly for e and £, e [0, Qo/2k], then one would obtain the desired estimate (by 
doing all the previous transformations backward): 

- I Xi\r-r*Heo]-^hs{f)fM-^dv>-C I X{\r-r*Heo]-Cwax{f)fM-^dv. 

To study the convolution integral of (13.91 1. we make a second splitting between \zu -v\ ^ 
g\ > and Irl; - i5| < (for some gi > 0). It gives 



\v-iv\Mi(w) 

dw 



a;eR2 {\V - a>|2 + ,^2)1/2 

r \v-w\Mm r Mm 

TB dw^ gi — dw 

Jm-v\>g,] {\V - W\2 + ef J{\w-v\>e,] {\V - Wf + 52)1/2 

,Jr MM r Mm ,A 

\ Jr2 {\v - ZDp + ^2)1/2 J{\^l,-zHg^] {\v - lI^P + 42)1/2 J 
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Then we use the obvious estimates 



Jk^ (Iz; - Z-y|2 + §2)1/2 1 + |p| 

for an expHcit constant C(k, Qq) > depending only on k, Qq, and 



3,Po/2k], I 



V i; e V ^ e [0, po/2k], I dzZ) < C(k, a) e"'"'' 



for an explicit constant C{k, Qi) > going to as gi goes to 0. It yields for small enough 
(depending on Qq) 

r Mi{w) ^ 1 r Mi{w) ^^^^^ 

dio^- dw. (3.10) 

for any v e ]R2^ § g [q^ Qo/2k] (we also refer to the Appendix A of this paper for a con- 
struction of the parameter gi). Consequently, for this choice of we obtain (|3.9l l with 
C = gi/2, i.e., 

- I X{\r-r-\<eo]£hs{f)fM~^dv> I X{\r-r*\<eo]-Cwax{f)fM~'^dv. 

Jr3 ^ Je? 

This, together with estimate ( I3.8l l, yield 

2)hs(/)^min{£,|}£»™x(/) 

which concludes the proof. □ 

Remark 3.6. The constant C* from the proof can be optimized according to the parameter 
go, by expliciting g\ as a function of go. Precisely, making use of Lemma A.l given in the 
Appendix, 

C = mm^ — , — > > 



.2k' 2 j ^5 

with rj = erf""*^ (^jj where erf~^ denotes the inverse error function, erf~^(|) 

0.4769. Notice that this lower bound for C* does not depend on the parameters a, |3. 

Remark 3.7. The above Proposition provides an estimate of the spectral gap of X.hs in 
L^{M~^). Precisely, we recall from \ \ \ that the spectrum of X,y\^ is made of continuous 
(essential) spectrum {A e IR; A < — Vq} where Vq = inf^gjjs 0hs(^) > and a decreasing 
sequence of real eigenvalues with finite algebraic multiplicities which unique possible clus- 
ter point is — Vq- Then, since is an eigenvalue of associated to M, one sees from the 
above Proposition that the spectral gap jUhs o/Xhs satisfies 

fZhs := min{A : -A e (-vo,0),-A e SCXhs) \ {0}} > C^max > '^""""'"'^^^'^ ~ ■ 

To summarize, one gets the following coercivity estimate for the Dirichlet form: 
Theorem 3.8. For Q = Xhs or Xmax, one has the following: 

- f Q{f){v)f{v)M-\v) dv > ^i\\f - gfM\\l,,M-^ "if e m) 

where, gf = | f{v) dv, and fi = fimax whenever Q = Xmax while, for hard-spheres 
interactions, i.e., Q = Xhs. one has > C*jUmax- 
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Proof. If Qf = 0, the proof follows directly from Proposition 13 . 3 1 and ( 13.71 ). Now, if / is a 
given function with non-zero mean Qf, seth - f - QfM. Then, p;, = so that 



■ Q{h)(v)h{v)M-\v) dv > 



This leads to the result since = <a(/) and <a(/)di; = 0. □ 

Adopting now the entropy language, one obtains the following relaxation rate, which is 
also new in the context of linear Boltzmann equation: 

Corollary 3.9. Let /o(c) e A1~^(c) di?) be given and let f{t) be the unique solution 

of dl.ll) with initial condition f{t = 0) = fo- Then, for any i ^ 0, one has the following 

- MIlhm-^) < exp(-f/0 ||/o - ' > 0, 

where [i = flmax when Q = Qmax while, for hard-spheres interactions, i.e., Q = (3hs. one 
has H > C*/imax- 



We state another corollary of the above Theorem l3.8l in which we strengthen the coerci- 
vity estimate: 

Corollary 3.10. For Q = Xtis or Xmax, there exists Cg > such that 

- f Q{f){v)f{v)M-Hv) dv > cMf - QfM) V/ e ^Q) 

where, Qf — I f{v) dv and a{v) is the collision frequency associated to Q. 
Jr3 

Proof. If Q - Xmax, since amax(^') = 1 the estimate is nothing but Theorem l3.8l Let us 
consider now the hard-spheres case, Q - Xhs- Arguing as in the proof of Theorem [33] it 
suffices to prove the result for f± M, i.e., whenever Qf = 0. We recall from 1 1 1 that Xhs 
splits as 

where J{f is a bounded (and compact) operator in L^{Ai~^). We then have 

\\f^^^s\\LM-.)= f -r{f)fM-'dv- f £Uf)fM-'dv 
Jr3 Jr3 

< W-rMfW^M-^. + %s{f) ^ + lWhs(/) 

\f'max^ / 

where \\J^\\ stands for the norm of as a bounded operator on L^{Ai~^) and we used 
Theorem l3.8l The corollary follows with 



11^11 + C>max' 



C,T = 



Remark 3.11. Here again, as in Prop. 15.51 the constant Cg > can be quantitatively 

estimated using for instance the estimate\\.J(^\\ < {h^r)^ 'sj^m^ ^^'^^ '^'^^ deduced without 
major difficulty from the explicit expression of J(f provided in [ 1] with t = (1 — 2k) Ik > 0. 

Remark 3.12. Recalling that ffhs behaves like (1 -I- \v\), the above corollary allows to 
control from below the entropy production functional by the weighted L^((l + |c|)Al~^, dv) 
norm. Such a weighted estimate shall be very useful for the diffusion approximation. 
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4. Diffusion Approximation. We shall assume again in this whole section that Mi = 0. 
From the results of the previous section, it is possible to derive some exact convergence 
results for the solution of the re-scaled linear kinetic Boltzmann equation 

edt Mt, x,v) + v V.,/,(f, X, v) = ^Qim, X, v), (4. 1) 

with initial condition /t.(x, i?, 0) = foix,v) ^ 0, with {x,v) e IR-' X IR''. Note that all the 
analysis we perform here is also valid if the spatial domain denotes the three-dimensional 
torus T^.. One shall prove that /f converges, as e ^ 0, to Ai{v)g where g = g{t,x) is the 
solution to the (parabolic) diffusion equation: 

dtg{t, x) = V:, ■ (DW^git, x) + uip), t > 0, x e E?, 

C (4 2) 

g{0,x) = goix) = fo{x,v)dv 

where the diffusion coefficient D depends on the model we investigate (hard-sphere inter- 
actions or Maxwell molecules). One shall adopt here the strategy of IJ2ul4J . Namely, to 
prove the convergence of the solution to ( 14.1b towards the solution g of (I4.2l l. the idea is to 
use the a priori estimate given by the production of entropy, as in f\A\ where this idea was 
applied to discrete velocity models of the Boltzmann equation. Let us define the number 
density and the current vector 

gi:{t,x)= \ fi:{t,x,v)dv, ic{t,x)=- f f^{t,x,v)v dv. 
Jr3 £ Jr? 



We also define he as 

h. 



e{t,x,v) = ^{fe{t,x,v)- ge(t,x)M{v)j. 



Integrating (14.11 ) with respect to x and v and using the fact that the mean of Q{fe) is zero, 
one gets the mass conservation identity 

I fe{x,v,t) dxdv - I fQ{x,v)dxdv, (4.3) 

Jr^xRJ JrJxR;! 

which means (using the fact that the equation preserves non-negativity) that, for any T > 0, 
the sequence g^ix, t) is bounded in L°°(0, T; L^(R^)). Now, multiplying ( 14. Il l by /t A4~^ and 
integrating over we get 

14: f fHt,x,v)M-\v)dxdv+ ^ f dwJvf,^{t,x,v))M-\v)dxdv 

-\ f feQ{fe)M-' dxdv = 0. (4.4) 

Now, because of the divergence form of the integrand, one sees that the second term in 
( I4.4l i is zero while, because of Corollarv l3.10l 

-\ f feQ{fe)M-^ dxdv 

^ Jk^xR?, 

- Ca I hf{t,x,v)M~^iv)a{v)dxdv. 

JrJxr?, 



INELASTIC SCATTERING BOLTZMANN MODELS 



17 



Consequently, Eq. (14.41) . together with (14.51) . leads to 



i-^ r f^(x,v,t)M~^{v)dxdv ^ -Ca [ h^^{x,v,t)a{v)M~'^{v)dxdv. 



2 

Defining therefore the following Hilbert space: 

= L^(ILl X E^,M~Hv) dxdv) 
endowed with its natural norm || • one has 

||/,(f)||' + 2c„ r \\h,is) V^lf ds < ||/o||' , Vi > 0. (4.6) 
Jo 

We obtain the following a priori bounds: 

Proposition 4.1. For any e > 0, let fi;{t) denotes the unique solution to ( 14.1b with fg e ^H, 
fo > 0. Then, for any < T < oo 

1. The sequence (fe)E is bounded in L°° (0, T; "H) , 

2. the sequence ( -sfahi^i; is bounded in (0, T; fH) , 

3. the density sequence is bounded in L°°(0, T ; L^(]R^) D L^(Rj)), 

4. the current sequence is bounded in [l'^{{0, T) X R^)] . 
Proof. The first two points are direct consequences of (14.61 ) with 

sup WMLioj-.-H) ^ \\Ml,v' sup II yi^h,\\ < {2c,r''^ ||/o|| . 

Now, Eq. ( 14.31 ) proves that the number density sequence {g£)c is bounded in L°°{0, T ; L^(]R^,)) 
and, according to Cauchy-Schwarz inequality, 

< Q,{t, X' v)M-\v) dv^ ' 

we see from point (1) that {gt)^ is also bounded in L°°(0, T; L^(]R^)). Finally, since = 
QcM + ehi; and J^, vAi{v) dv = 0, one has 

Jdt I |;f(t,x)|^ dx = I df I dx I vhi;{t,x,v)dv 
Jul Jo JrI JrI 

while, from Cauchy-Schwarz inequality and the fact that a is bounded from below 

If vh,it,x,v)dv f \v\^M{v)dv]( f hfM-^dv] 

so that ^ ^ 

[ dt f \i,{t,xf dx^— f WhMllvdt 
Jo Jv.y m Jo 

and the conclusion follows from point (2). □ 

Remark 4.2. Since = ehi, + g^M, noticing that ^^r^ a{v)A\{v) dv < oo, one deduces 
from the above points (2) and (3) and that the sequence ( 's/ofi)^ is bounded in (0, T; "H) . 

For any T > 0, we define 

Qr = (0, T)xUlx and d^j = dx dv dt. 
The bounds provided by Prop. l4.ll allows to assume that, up to a subsequence, 

/, ^/ in L^{Qt ; oM'^ dfir), h, ^ h in L^inr ; oM'^ dfij); 
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g,^g in 1^(0,7) xR% j, j in [l\{0,T)xRI)^ . 

Let^ e L^{Q.T, ff~^Aldfir) = [L^(Qr/ aM~'^ dfij)] be given. Since a = ffmax is constant 
while a - dhs behaves asymptotically like (1 + one easily has from Cauchy-Schwarz 

<p{t,x)= [ M{v)^{t,x,v)dveL'^{{0,T)xK^), 

and therefore 

lim I dt I gi:{t,x)(p{t,x)dx = I dt I g{t,x)(p{t,x)dx. 
Thus, writing /t = f(f Al + ehe, one checks that 

lim I /£^d/,ir = I di I g{t,x)(p{t,x) dx = I fiAIWdfir, 
i.e., /c ^ gAi in L^(i2r, ctA1~^ dfir)- In particular, f{t, x, v) = g{t, x)M{v). Moreover, 

lim f h,W dfiT = I hW dfir- (4.7) 



for any W = W{t, x, v) e L^{Qt, o d/.ir). Now, choosing W independent of v, one sees 
that 

I h{t, X, v) dv = 0, Vi > 0, X e R^,. 

Finally, using in ( 14.71 ) a test function i^) = v(p{t,x) with (p e L^((0, T) X ]R^), we 

deduces from the weak convergence of to that 

i{t, x) - I vh{t, X, v) dv. 

Finally, integrating equation ( 14.11 ) over leads to the continuity equation 

dtg,{t, ^) + div:,;;(f, x) = 0, Ve > 0. (4.8) 

We deduce at the limit that 

dtg{t, x) + divxjit, x) = 0, t > 0, xeT^ (4.9) 

in the distributional sense. We summarize these first results in the following: 

Proposition 4.3. Under the assumptions of Proposition \4~l\ for any T > 0, up to a subse- 
quence, 

i) (gc) converges weakly in L^((0, T) X ]R^,) to some g; 

ii) (h^) converges weakly in L?{^1t, a Ai~^ dfj-j) to some function h with 

I h{t, X, v) dv = 0; 

Hi) (fe) converges weakly to gAi in L^{Qt/ aM~^ d\ij); 

iv) (jc) converges weakly to i{t,x) = \ vh{t,x,v)dv in\L^{{0,T)xE.l)\ . 

where g and j are related by i4.9\ . 

The problem of the diffusion approximation is then reduced to the one of finding a 
suitable relation, similar to the classical Pick's law, linking the current j{t, x) to the gradient 
of the density g{t,x). Such a Fick's law (and the corresponding coefficient) shall depend 
heavily on the collision kernel. 
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4.1. Maxwell molecules. When dealing with Maxwell molecules, i.e., whenever Q = 
Xmax, it is possible to obtain an explicit expression for the diffusion coefficient. Precisely, 
multiplying equation ( 14. Il l by v and integrating over gives 

£^dti,{t,x) + l f {v^v):y,Mt,x,v)dv]=- [ £maAfe)vdv (4.10) 

Now, as we already saw it (see (13. 2t ): 

■£max{fe)'vdv = -a(l - j3)e;, = -Ao,ie;,. 



X?, 



Then, recalling that /^ (f, x, v) = g^{t, x)M{v) + ehdt, x, v), Eq. ( 14.101 ) becomes 

£^dtie{t,x) + A:V,g,{t,x) + £i i (v®v) :V:,h,{t,x,v)dv] = -Ao,i je (4.11) 



where A is the matrix of directional temperatures associated to the distribution Ai: 
A = I {v®v) M{v) dv - — Id - diag 



m \ m m m 

One may rewrite ( 14.1 11 1 as 
m 



-V^g,{t,x) + Ao,iii:{t,x) = -£^dti,{t,x)- e^J ^{vigiv) : V:,h,{t,x,v)dvj 



Choosing a test-function ijj e '^^{{0, T) X R^), the above equation reads in its distributional 
form: 



0# 
m 



f dt f V ^ip{t, x)g,{t, x) dx - Ao,i f dt f 4>{t,x)i,{t,x)dt = 
Jo JkJ ' Jo JrJ 

- I dt \ dtip{t,x)ji:{t,x)dx- £ I Z7) : V.vi/'(f,x)dfi 

Jo Jr? Jor 



T 

and, by virtue of the bounds in Prop. 14. ll the right-hand side converges to zero as e ^ 
and one gets at the limit: 

i{t,x) = —^V,g{t,x) (4.12) 

in the distributional sense. The above formula provides the so-called Pick's law for Maxwell' s 
molecules. One deduces the following Theorem: 

Theorem 4.4. Let foefi and, for any £ > 0, let f^t, x, v) denotes the solution to ( 14. 11 1. 
Then, for any T > 0, up to a sequence, converges strongly in L^^^{ilT ; Ai~^ d^fr) to- 
wards g{t, x)M{v), where g{t, x) is the solution to the diffusion equation 

d,g = y,-\-^Vxg{t,x)], g{t = 0,x)= f Mx,v)dv. (4.13) 
\mAo,i I J]R3 

Proof. We already proved that converges weakly to gAi in L^((0, T) ; 'H). To prove the 
strong convergence, since 

r \\Mt)-Q.{t,x)M\\l, = £^ [ WhMl'n^O 
Jo Jo 

it suffices to prove that g^Ai converges strongly to gAi in 'H. This is equivalent to prove 

that converges strongly to g in L^(0, T ; L^^^(]R^)). This is done in the spirit of L14| and 

lfT2l by using a compensated-compactness argument. Precisely, let us define the following 
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vectors of xM^ : Uf = (/t, Q^) and Vf = (0, g^). From ( 14.81 ), one sees that diVj:,fU£ = 0, 
in particular (diVx,tUt)t is bounded in L^(]R^, X ]R^). Now, from ( 14.111 ). one sees that 
A : VvOf IS a bounded family in L2((0, T) X M^). 

Since A = —Id, it is clear that 

curlV -I -'^■^•^^ 
' \ V.,^), 

is bounded in [Lf„^((0, T) X IR^)]*""^. Now, from the div-curl Lemma ElEa, U, • = g^^ 
converges to in D'{{0, T) X R^). 

Moreover, we already saw that g, is bounded in L°°(0, T ; L2(]R3)) from which we deduce 
the strong convergence of g^ to g in L^((0, T) ; L^^^(]R^)). □ 

Remark 4.5. As already pointed out in |211, f/ie dependence of the dijfusivity Dmax 
0*/mAo,i on the inelasticity parameter j} shows that inelasticity tends to slow down the 
diffusive process. 

4.2. Hard spheres. When dealing with hard-spheres interactions, it appears difficult to 
obtain an explicit expression of the diffusion coefficient. Nevertheless, its existence can be 
deduced from Theorem [33] Indeed, a direct consequence of the Fredholm Alternative is 
the following: 

Proposition 4.6. For any i = 1,2,3, the equation 

£hs{xd = ViM{v), veE? 

has a unique solution Xi G L'^{a{v)A\~^{v) dv), such that {Xi,M) = I Xi{'^)d'v = Ofor 

M 

any i = 1,2,3. 



Remark 4.7. Note that the above Proposition holds true only because we assumed the bulk 
velocity iii to be zero, i.e., J^3 vAXAv = 0. If one deals with a non-zero bulk velocity U\, 

then if one denotes a{v) = v — Ui, Xi then solves Xhste) = ai{v)M{v) (see also ( 14.151 )), and 
moreover the limit diffusion equation includes in this case an additional drift term Ui ■ V^p, 
see Eq. ( l4!2b . 



Then, setting x = (Xii Xir Xi) one defines the diffusion matrix: 
D := - I v®x{'v)dveU^''^. 



Adapting the result of f20l, the diffusion matrix is given by D = diag(Dhs, Dhs, Dhs) for 
some positive constant Dhs > 0, namely. 



Dhs = - ( vixi{y)dv = - I -Chs(XikiA1 ^ dz; ; 

Jv.f. JtrI 



l^Wx 



||2 

iIl2(m-i)- 



Remark 4.8. Note that, when dealing with Maxwell molecules, for any i = 1,2,3, the 
function Xi appearing in Proposition \4.6\ is given by Xi — and we find again the 



expression of the diffusion matrix D = ^j^— Id. 



Recall that, for any T > 0, we defined Or = (0, r)x]R,3x]R3 and d^ij = dx dv dt. Then, 
for any <p e L°°(]R^) and any ^p e '^^{{0, T) x W^l), multiplying Eq. (|4l]) by (p{v)ijj{t,x) 



INELASTIC SCATTERING BOLTZMANN MODELS 



21 



and integrating over Qr one has 

I g,(t,x)M{v){v ■y^\p{t,x))(p{v)diUT + I £hs{K)cj){v)4>{t,x)di.iT = 

e{ f (Pf,dt4> d^iT+ [ he{v-V,ip)(pdfiT\ (4.14) 



In particular, by virtue of Propositions RTTI and K. 31 one sees that 



lim ( I M{v)g,{t, x) (v ■ Vxip{t, x)) (p{v) djUj + \ £hsihE)(p{v)ip{t, x) d^r] = 0. 
f ^0 V JsIt Jht > 

)w, one deduces easily as in [121 that 

r Q{t,x){vV:cW>^))M{v)(p{v)d^T = - \ £t,s(h)(p{v)ip{t, x) djlT, (4.15) 



which means that, in the distributional sense, 

div:,{vM{v)g{t, x)) = Xhs(^), f > 0, X e R^. 
Since h is of zero R^-average, Proposition l4.6l asserts that 

h{t, X, v) = -xiv) ■ yxQ{t, x) 
and Proposition |4.3l (/v) leads to 



,x) = I vh{t, x,v)dv = Y) : VxQ{t, x). 
Jr3 



We then obtain the following: 

Theorem 4.9. Let < fo{x, v) e L^(]R^xR^, dv) be given and let fs, be the associated 
sequence of solution to ( 14. Il l where Q = Qhs- Then, up to a subsequence, converges 
strongly in L^^J^Clj,M~^ dfir) to g{t,x)A\ where g is the solution to the parabolic 
diffusion equation (14.21) where the diffusion coefficient Dhs is given by 

Dhs-- r 1^1X1(1^) dz;e]R3x3 

with Xi defined in Prop. \4.6\ 



Proof. We aheady proved that /t converges weakly to gAi in L^((0, T), 'H) and the strategy 
to prove the strong convergence is that used in Theorem 14.41 Precisely, we define again 
Uf = (/f , ) and Vt = (0, g^) and observes that again (diVv,fUf )e is bounded in L^(]R^ X 

Rj^). Now, from (I4.14l i. with (p{v) = — and setting T = I M{v) dv, one sees that 

g^ satisfies: 

r:V.^^)f= £^s{hi)—dv - e\dt I — dz; + div.,. ——h^dv 

so that r : Vj-ff lies in a bounded subset of lfo,((0, T) x R^). Since F is invertible, one 
proceeds as in the proof of Theorem 14.41 that g^ converges strongly to g in L^^^((0, T) X 
R3). " □ 

As we saw it, the diffusivity D^s associated to hard-spheres interactions is not explicitly 
computable, the solution x not being explicit. It is however possible to obtain a quantitative 
estimate of Dhs in terms of known quantities {i.e., that do not involve Xi)'- 
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Proposition 4.10. One has the following estimate: 



< Dhs < 



Chsm Ao,iC*m 



where C* is the constant provided by Prop, \3.3\ and Chs = ^ > 0. 

Proof. We begin with the lower bound of Dhs- For any s e R, let 
f (s) = {£hs{Xi + sviMlxi + sviM). 
Since Xhs is negative, one has P{s) < for any s e ]R. Moreover, 

P{s) = (Xhs(Xi),^i) + 2s<Xhs(Xi),i'iM) + s^{£t,,iviM),v^M) 
= -Dhs + 2s||i;iM||^,(^_i) + s\£us{viM),ViM). 
We get therefore that 

Dhs > 2s\\viM\\%^j^_,^ + s^UusiviM), viM), Vs e R. 
With the definition of Chs (note that Chs > since Xhs is negative and ViAi±M), we get 

Dhs > (2s - Chss^) \\viM\\l,^j^_,y Vs e R. 
Optimizing with respect to s, one sees that 

Dhs > —WviMW'. = 

Chs ^ ' ChsJW 

To get an upper bound for Dhs, we use the fact that, thanks to (I3.2l i. 

Dhs = -{Xl,1^lM) = Aol{£vnax{Xl),^lM). 

Now, as above, for any s e R, define Q(s) = {£max{sxi + ^iM),sxi + ViM). Here again, 
Q(s) < for any s e R and 

Q(s) = S^{£max{X\),X\) +2s{£^ax{X\),V\M) + {£max{viM),ViM) 
= s^{£vnax{Xi),Xi) + 2Ao,iDhsS - Ao,i||t7iAl||^2(^_i). 
Now, according to Prop. |331 (XmaxCxO/A:!) > ^{£hs{Xi)> Xi) = -Dhs/C* so that 

> Q(s) > -^s2 + 2Ao,iDhsS - Ao,i\\viM\\l,^j^_,y Vs e R. 
Optimizing the right-hand side with respect to s e R, we get 

Ao,iII^iM|Il2(m-i) > A^^C^Dhs 
which gives the desired upper bound. □ 

Remark 4.11. It is possible to provide some upper bound for Chs- Namely, using the fact 
that there exists Vi > such that a{v) < Vi(l + li^l), it is easy to see that 

'^^^^n^i r o{vfv\M{v)dv^^ r {l + \v\fM{v)dv. 

This very rough estimate could certainly be strengthened. Note also that the upper bound 
for Dhs reads as 

t(i - a)(i - |S) v^;^^ 



Dhs < 



a(l - jS)(l - a(l - j3)) m 
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V5 

where t = — 3.3154 is a numerical constant and we used the lower bound 

erri(l/2)V2 
q/C* provided by Remark \3.6\ 



Appendix A. We provide here a constructive proof of the coefficient Qi appearing in 
Proposition 13.31 with the aim of finding quantitative estimates for the coefficient C* in 
Prop. 13.31 Namely, recalling that k - a{l - |S), one has 

Lemma A.l. Given qo > 0, there exists gi > such that 

r Mm ^ . ^ 1 f Mm 

Jiiffi-f.i<pi) {\v - + B-r 2 -J^^ (1^ - + i^r 

for any v e IR^, (5 e [0, ^o/2k]. Moreover setting rj = ^2§Lerf~^(i), where erf~^ is the 
inverse error function, one has 



4k2 



VO < po < 2kj]. 



Proof. We assume without loss of generality that Mj = 0. Let > be given. Let us fix 
V - (01, V2) € and £, e [0, ^o/2k]. Using polar coordinates, it is clear that 



\2n0j J|ur,_„| 



Mm 

aw 



(|0 - II'P + 52)^^^ 

, fft rexp(-fl/-2) f2" . . 
exp(-fl|z?r) — dr exp - 2flr(z;i cos + sin 0) d0 

Jo V^^n2 Jo ^ ^ 

where a = »ii/(20i). Therefore, a sufficient condition (independent of v) for (lA.ll i to hold 
is that 

r^' rexp(-flr2) l r°° rexp(-flr2) 

— m^dr<- Zl— ^ dn V5e[0,^.o/2K]. 

-'0 + ^2 Z Jo ^^-2 _,_ ^2 

It is not difficult to see that this is equivalent to 

eri(^^a{gl + e)) - erf( V^5) < ^ - ^erf( V^5), e [0, Qo/Ik] 

where erf is the error function erf(x) = J" exp(-t^) dt, x > 0. This allows to define a 
function: 

z : 4 e ]R+ 1-^ z{£.) 
where z(^) is the nonnegative solution to the identity 

V«(z' + = erf-i + ^erf( V^5)) (A.2) 

where erf~^ is the inverse error function. Clearly the Lemma is proven provided 
pi ■- min{z(.S), £. e [0,go/2K]] > 0. 



Note that, according to ( IA.2I ). the function z( ) is continuously differentiable and there is 
some C e [0, f)o/2K] such that mm{z{£.), e [0, po/2k]} = z{Q. In particular z'(C) = 
and one checks, thanks to ( IA.2b . that 



1 

z'{E,) = - Vz2 + e exp(flz2) - 4, V£ ^ 0. 
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In particular, z'{Q = is equivalent to 

4C'exp(-flz2(Q) = z2(Q + c2, (A.3) 

and one sees that z^{Q = should imply C = whereas, according to ( IA.2I ). z(0) 0. 
Consequently, all the local extrema of z are positive. Therefore, 

gi = z(C) = min{z(<5), < ^ < go /2k] > 

which achieves to prove that (lA.ll) holds true for some gi > 0. It remains now to provide 
some estimate for gi . Precisely, defining 

we see from iA.2\ that z^{^) + ^ if, for any (J ^ 0, so that ^ rf - for any 
^0 G (0, 2k/]), which achieves to prove the lemma. □ 

Remark A.l. According to the above Lemma, with the choice of go = one obtains 
r/.armin(f,f)>^. 
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